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Abstract 

The shape of the neutrino pulse from the supernova SN1987a provides one 
of the most stringent constraints on the size of large, compact, "gravity- 
only" extra dimensions. Previously, calculations have been carried out for 
a newly-born proto-neutron star with a temperature of about 50 MeV at 
nuclear matter density. It is arguable that, due to the extreme conditions 
in the interior of the star, matter might be a quark-gluon plasma, where the 
relevant degrees of freedom are quarks and gluons rather than nucleons. We 
consider an energy-loss scenario where seconds after rebounce the core of the 
star consists of a hot and dense quark-gluon plasma. Adopting a simplified 
model of the plasma we derive the necessary energy-loss formulae in the soft- 
radiation limit. The emissivity is found to be comparable to the one for 
nuclear matter and bounds on the radius of extra dimensions are similar to 
those found previously from nuclear matter calculations. 
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I. INTRODUCTION 



Ideas that we live in a world consisting of d compact extra dimensions in addition to 
the usual four infinite dimensions are not new and date back to the 1930s. They imply 
that there are Kaluza-Klein (KK)-modes corresponding to excitations of ordinary standard 
model particles in the extra dimensions. Such modes, however, have not been seen in any 
collider experiment so far Q. 

Recently, a variation of this concept has been revitalized by Arkani-Hamed et al. H-fl] 
who considered an alternative picture in which standard model fields are confined to a 
four- dimensional "brane" while gravity, constituting the dynamics of space-time itself, is 
allowed to propagate into the whole (4+<i) dimensional "bulk". This picture caused some 
excitement because it provides a natural solution to the hierarchy problem. As detailed 
in ||, the extension of the four-dimensional space for gravitons leads to a "dilution" of 
gravity into the extra dimensions and therefore to the smallness of Newton's gravitational 
constant Gn- Matching the Planck mass Mpu^ in the four-dimensional world to the one in 
the (4+o?) dimensional world, M P i( 4+(i ), via the relation 

Ml m ~ R d M^ +d) , (1.1) 

one can bring the Planck mass down from 10 19 TeV to the weak scale of 1 TeV by choosing 
an appropriate size R of the extra dimensions.^] For one extra dimension, R ~ 10 11 m, in 
this case gravity — specifically Newton's inverse-square law — would be modified on the scale 
of astronomical distances, a case clearly ruled out. But already for d = 2, R ~ 1 mm, which 
is just out of reach to present experimental measurements of the gravitational force law ||. 

Another way of testing these models is to search for "beyond the standard model" physics 
at the scale M P1 ( 4+rf ). Even though one would naively expect a violation of the standard 
model at this scale to be of only gravitational strength, the large number of excited KK- 
modes at such high energies compensates for the weak coupling between these KK-modes 
and the ordinary matter fields. Such effects could alter standard model predictions due 
to virtual KK-mode contributions or manifest themselves in the form of massive graviton 
production 10-0. It might be possible that such effects can be seen in the form of missing 
energy at future colliders operating at the TeV scale like the CERN LHC ||10|.[TT[ . 



The most stringent constraints on the size of extra dimensions come from astrophysical 
considerations. Recently it has been pointed out that an overproduction of massive KK 
modes in the early universe might result in early matter domination and therefore a lower 
age of the universe |12 . 



Another way to obtain bounds on the size of extra dimensions from astrophysics comes 
from supernova SN1987a. Our current theory of supernovae predicts the shape and duration 
of SN1987a's neutrino pulse very well. Hence, if some new channel transports too much en- 
ergy from the interior of the supernova then the current understanding of SN1987a's neutrino 



1 A different way of solving the hierarchy problem has been proposed by Randall and Sundrum ||. 
In their model the metric is non-factorizable but rather the four-dimensional metric is multiplied 
by an exponential "warp" factor. We do not consider this model here. 
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signal gets invalidated. Raffelt |L3] has quantified the maximum possible emissivity for a new 
energy loss process which doesn't conflict with our current understanding of the neutrino 
signal as 10 19 ergs/g/s. Using this criteria there have been several calculations [|T4]-|i6| to ob- 
tain bounds on the size of the extra dimensions for a newly born proto-neutron star assuming 
that the interior of the star consists of nuclear matter at a temperature of ~ 50 MeV. A 



more rigorous approach has been undertaken recently by Hanhart et al. |17| who performed 
detailed simulations of the effect of exotic radiation on the neutrino signal. 

Although the star's inner core might well consist of nuclear matter, because of the extreme 
conditions seconds after core bounce with a density in some regions possibly up to ~ 15 times 
the nuclear matter density n = 0.17 fm -3 and a temperature of maybe up to ~ 100 MeV, 
matter is near a phase transition to a deconfmed QCD plasma. 

It is therefore important to investigate how bounds on the size of extra dimensions change 
if a scenario is considered where, in contrast to taking nucleons as the relevant degrees 
of freedom, the calculation is carried out using quarks and gluons as degrees of freedom. 
However, in the regime we are considering here matter is not in a state where the density 
and temperature are so high that a first order perturbative calculation would suffice; it is 
rather in a regime where deconfinement just sets in and the QCD coupling constant is still 
large (a s ~ 1). For this reason we will carry out our calculations for a (unphysical) regime of 
very high density (up to 10 7 n ). There, the coupling is weak and perturbative calculations 
are feasible. From there we extrapolate down to get the emissivity in the physical region. 

Note that if the temperature falls below a critical temperature T c of ~ 50 MeV (for 
a density of a few times nuclear density) this very dense plasma most likely undergoes a 
phase-transition to a color superconducting phase (for a recent review see [|18j). Since we 
do not consider such a scenario here, we assume the temperature to be above T c . 

The paper is organized as follows: In Section [TT] we calculate the amplitude for the 
emission of soft gravi- and dilastrahlung into KK-modes from quark-quark (qq) scattering 
in a degenerate QCD plasma. This calculation is fairly general, it applies to soft radiation of 
KK-modes from light relativistic fermions. Then, in Section |TJ, we consider a simple model 
of the plasma consisting of three light quark flavors to estimate the size of the qq scattering 
cross sections. To treat the occurring divergences for colinear scattering we incorporate the 
effects of the surrounding plasma into the gluon propagator by introducing a cutoff mass. 
We do not consider any other many-body effects; our goal is to obtain an estimate for the 



qq scattering cross section, an exact calculation is presently not feasible. In Section [TV 



we carry out the phase-space integration to get a formula for numerical calculation of the 
emissivity of a gas of weakly interacting quarks. We also provide an approximation formula 
for the degenerate limit. Finally, in Section M, we calculate the emissivity and the bounds 
on the size of the d extra dimensions for the case of d = 2 and d = 3 and compare them to 
the result from previous nuclear matter calculations. 



II. 4D-GRAVITON, KK-GRAVITON, AND KK-DILATON RADIATION FROM 

QUARK-QUARK SCATTERING 

We assume the d extra dimensions to form a <i-torus with radius R so that the KK- 
graviton and KK-dilation mode expansions are given by 
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and (p ab (x, y) = J2 h ab( x ) ex P * 



j ■ y 

R . 



(2.1) 



respectively. In addition to these particles there emerge massive spin-1 particles from the 
KK reduction of the Fierz-Pauli Lagrangian which decouple from matter and need not 
be considered here.0 

The coupling of the 4D-graviton g = h, KK-graviton h\ and KK-dilaton J to the 
energy-momentum tensor T^ v of the matter field is governed by the Lagrange density 



3{d + 2y aa " 



(2.2) 



where k is defined as k = vI&ttGjv with Gn being Newton's gravitational constant. Note 
that in Eq. fl2.2|) the massless graviton — being the zero-mode of the modes — has 
been written separately for clarity. The energy-momentum tensor T^ v for a free quark is 
given by 



T, 



(2.3) 



As we will discuss in the next section, the leading contribution from gravitational emission 
for the degenerate quark-gluon plasma of our model comes from the qq scattering processes 
qq — > qqX, where a low-momentum 4D-graviton (X = g), KK-graviton (X = h), or KK- 
dilaton (X = 0) is emitted. 



A. Amplitude for qq — > qqX in the Soft Limit 

For a degenerate plasma the participating quarks are near the Fermi surface and hence 
the momentum k of the emitted graviton is small compared to the quark momenta. In 
the soft limit the leading diagrams for the process qq —>■ qqX are those where the graviton 
radiates off one of the external legs as shown in Fig. |l] (bremsstrahlung process). For small 
k these diagrams are 0{k~ 1 ) while diagrams where the graviton couples to some internal 
line are 0(1) and therefore suppressed |2T)|-[2"2" 



The amplitude M^ v of this process then factorizes into the amplitude A (which includes 
a sum over spins and colors of the incoming and outgoing quarks) for on-shell qq — > qq 
scattering and the conserved gravitational current T^ v . To lowest order in k, M^ v can be 
written as [11(3] 



Kt l 1 



^=«E m^—- A, (2.4) 
z j=i Pi ' ^ 



2 Note that while the torodial compactification is conceptually simple, it might not be realistic 
since it doesn't take into account curvature ("warping") caused by fields in the bulk and on the 
brane. As has been recently investigated by Fox [O], weak warping might increase the emissivity 
by a factor ~ 2 and therefore strengthen the bounds. 
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FIG. 1. The leading order diagrams for the process qq — > qqX. In the soft limit these dia- 
grams are dominant over those where the graviton couples to some internal line. A is the on-shell 
amplitude for quark-quark scattering. 



where 




1,2 
3,4 



(2.5) 



Here, the pj denote the four- momenta of the incoming (j = 1,2) and outgoing (j = 3,4) 
quarks while k — (a?, k) is the graviton's four- momentum. The result Eq. Q2.4|) has been cal- 
culated by Weinberg |23[] using an entirely classical derivation. This is not unexpected since 
to lowest order in k graphs where the radiation comes from an internal line are suppressed 
and so the radiation from the external legs is the dominant one, corresponding to classical 
bremsstrahlung from the gravitational charges accelerated during the collision. 



B. Energy-Loss Formulae 

The energy loss into a KK-mode j and into the momentum space volume element d 3 k 
due to 4D-graviton, KK-graviton, and KK-dilation emission is given by 



d 3 k 



£ [M a ?e a p(k, A)]* Af"V(£, A) = to 



d 3 k 



(2vr) 3 2w 



M Xv *M Xu 



M 



(2.6) 



5 r 



d 3 k 

dei = uo- — — — > 

h (27r) 3 2c; ^ 



A=l 



(2.7) 



and 

dl 



d 3 k 



n(ra-l)/2 f 



10- 



3(d + 2) (2tt) 3 2cj 



A=l 



4(M) 



eUk,X) 



d 3 k 



■M„ a M, 



2{d-l) 



2(2vr) 3 q " 3(d+2)' 

(2.8) 
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respectively. Here, e^, e^,, and el v are the polarization tensors of the 4D-graviton, KK- 
graviton, and KK-dilaton, respectively 0. 

Utilizing the conservation of the energy-momentum tensor, k^M^ = 0, M^ u can be 
expressed in terms of the purely spacelike components as 

y y y 

M 0i = Mji, M 00 = -rMj-j (2.9) 

or, upon defining 

W^ij = ■q ilx r] jv - 77^7^ — - i] ^r] ju — + r] 0fl r] 0u ^-, (2.10) 

be written as 

M, u = W^M'K (2.11) 



Using this notation the energy loss formulae in Eqs. ( |2.6| )-( |2TE| ) can be written solely in 
terms of the space-space components of the tensor M^ v as 



di = ^S M « A$WMw ' (2 - 12) 



where the A 1 ^ 1 are given by 

Af = 5 ik 5 jl - -5 ij 5 kl - — Uk j k l 5 ik - k k k l 5 ij - k i k j 5 kl ) + -^-k' l k j k k k l , (2.13) 

11 9 c ).k 2 ' V / 9k -4 



AT = W^B^pW 



a(3kl 



S ik S jl + S il S jk _ _ 5 ij 5 kl} + k i k j k k k 



( <f fc FA;' + <P"W + 5 il Vk k + 6 jl k i k k - -(k l k J S kl + k k k l 5 ir ) ) , (2.14) 



and 



aT 1 = w aij w R /3kl ^^ — — 

* a * 3(d + 2) 



+ 



+ ±-k*yk k k 1 ^ 2{ - d 1 

CJ 4 



3(d + 2) 
(2.15) 



The tensor Bn Ua /3, defined by the spin sum of the KK-graviton polarization tensors 



A=l 



(2.16) 



is given in |7|f|. 

We will now choose the center-of-momentum (COM) frame as our reference frame where 
the 4-momenta of the two incoming quarks are given by (p°,p) and (p°, —p) and those of 



6 



the outgoing quarks by (p'°,p') and (p'°, —p'). The COM scattering angle 9 cm is defined by 
p. p' — pp' cos6 cra . We also assume the quarks to have equal mass. 

If the momenta of the scattering particles are much smaller than their rest mass M — as 
is the case in soft radiation from nucleon-nucleon scattering — it can be neglected in the 
denominator in Eq. ( |2.4|) and is simply given by |16 



M ij = 1 ^j( Pi p j -p' i p' j ) A. (2.17) 

In our case, however, this approximation is inappropriate since the quark mass is very small 
compared to the quark Fermi momenta, we therefore write Eq. ( |2.4|) as 




PiPj PiP'j PiP'j 



p°u + p ■ k p'°uj — p' ■ k p'°uj + p' ■ k , 



A. (2.18) 



The / dk/Air integration in Eq. ( |2.12j) can then be carried out analytically for the three cases. 
For soft graviton radiation the quark momenta are assumed to lie on the Fermi surface and 
it is justified to assume \p\ ~ \p'\. We can then approximate the momenta p and p' by 

f = P ^f. (2.19) 

Replacing p and p' with p the formulae simplify considerably and the / dk/47c integration 
can be written as 

/ ^M ij Kf l M kl = ^g x (x, y, 9 cm ) \A\ 2 (2.20) 
J Att uok 

in terms of A and the dimensionless functions gx which are given in the Appendix. Here 
we have also defined x = mj/u and y = p°/p = ^p 1 + M 2 /p. 

Note that for massless quarks one of the denominators in Eq. ( |2.18| ) becomes zero if 
the momentum of that quark is parallel to k and one might think that this term becomes 
infinite. However, as already noted in Ref. |23[ for the case of the massless 4D-graviton, this 



singularity is only spurious and gets canceled by A x kl in the numerator so that gx is finite 
not only for X = g, but also for X = h and X = 0. 

We can then employ the energy-momentum relation for the KK-mode j and write the 
formula for the energy loss per unit frequency interval into the given KK-mode of mass 
vr?-, = u 2 — k 2 as 



3 

J ,.,U r At _ _ -2^2 



de^ = uk_ fdk^ Kijia ^ _ K 2 p 
du 4tc 2 



J ^-M t] Af l M kl = ^JL gx ( Xiyi e cm ) \A\ 2 . (2.21) 



This result is general as it applies to soft radiation from any light relativistic fermions with 
equal mass, not just to the quarks we are about to consider. 
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III. A SIMPLE PLASMA MODEL 



After bounce the homologous core of the star undergoing a supernova explosion gets 
compressed to a very high density. Previous calculations [|14|-|T6|1 assumed a density of 
1 — 3 times nuclear matter density n = 0.17 fm and a temperature of T « 30 — 80 MeV 
corresponding to a (non-relativistic) nucleon chemical potential of [in ~ 50 — 100 MeV. 

Under such extreme conditions the star's core not only consists of free neutrons and 
protons but these degrees of freedom are highly excited and near a phase transition to a 
deconfmed quark-gluon plasma. It is therefore worthwhile to investigate how the graviton 
emissivity and hence the bounds on extra dimensions change if one carries out the calcu- 
lations in this regime. Our intention here is to investigate graviton radiation into extra 
dimension from such a QCD plasma. 

The physics of a QCD plasma under such violent conditions is particularly rich. In addi- 
tion to quarks of different flavors, a variety of other particles are present and being produced 
constantly, such as electrons, positrons, neutrinos, and collective excitations like plasmons. 
The properties of the quark-gluon plasma are far from being completely understood so that 
we do not even attempt to take the whole richness of matter under such extreme conditions 
into account here. Instead, we assume a rather simply model for the plasma (not unlike the 
one used in Ref. where, in addition to electrons and neutrinos, only u, d, and s quarks 



are assumed to be present. The possible weak reactions in the plasma are then 

d <-> u + e~ + v~ and s <-> u + + v~ . (3.1) 

At thermal equilibrium the chemical potentials have to add up to zero. For a low temperature 
T this requirement reads 

= V-u + + and fi s = n u + fi e + fi v . (3.2) 

In addition, we want to satisfy charge neutrality, thereby imposing the constraint 

2 11 

= Q = ~n u - -n d - -n s - n e . (3.3) 

These three conditions can be simultaneously satisfied by choosing 

n u = n d = n s = y = n B and n e = n u = 0, (3.4) 

where n q is the quark number density and % denotes the baryon number density. We do 
not consider dynamical production of these or any heavier quark flavors. As an additional 
simplification we take the mass of the u, d, and s quarks to be zero. 

The main contribution to graviton radiation from this plasma comes from the 
bremsstrahlung reaction qq —>■ qqX (X = g, h, <fi) where a graviton is radiated from an 
external quark leg. This reaction is of order 0(Ka s ) where a s is the QCD coupling con- 
stant. As we will see below, a s is rather large and 0(1). Moreover, we get contributions 
from processes involving dynamical produced electrons (e) and positrons (p), like ee — > eeX 
or ep epX. These reactions are suppressed compared to qq — > qqX by a factor of 
aQED/^s, where «qed = 1/137 is the fine-structure constant. More contributions come 
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from graviton-production reactions like q — > qXX. These, however, are 0(k 2 ) and therefore 
highly suppressed. We also assume that because of T <C \i q the number density of gluons is 
much smaller than that of quarks and therefore we do not consider the reaction qg — > qX. 
Note that the inclusion of any KK-modes generated by channels we neglect will only increase 
the emissivity and strengthen the bounds obtained below. 

The baryon number density ng in the central core region of a stable neutron star is 
~ 1 — 15 times the nuclear matter density no = 0.17 fm -3 , depending on the nuclear matter 
equation of state used (compare Ref. |24| and values cited therein). Assuming quark matter, 
the (relativistic) quark chemical potential fi q can be related to the quark density n q using 
the formula 



^ = 6 ^/ dee 2 N q (e) 

Z7T Z JO 



(3.5) 



where N q (e) = (exp[(e — fiq)/T] + 1) _1 is the mean occupation number and e is the (rel- 
ativistic) quark energy. Choosing n q and T we can then determine \i q \ for T = 50 MeV 
and n q = 10 — 20uq we find fi q ~ 500 — 600 MeV, which is in good agreement with other 
calculations [2^|. Therefore quark matter has /x 9 /T ^> 1 which confirms our assumption 
that it is degenerate. 

Furthermore, using fi q and T we can calculate the average momentum square (p 2 ) carried 
by the quarks using 



(P 2 ) 



Jd 3 pp 2 N q {e) 
Jd 3 pN q (e) 



(3.6) 



For typical parameters for the plasma (n q ~ 10 — 20no, T ~ 50 MeV) we calculate the 
average momentum carried by the quarks to be p pa 420 — 520 MeV which is just about 
the order of energy at which the strong interaction becomes strong and the QCD coupling 
constant a s is expected to be rather large. Of course, this is eventually a sign that matter 
under these conditions is at the phase transition to the QCD phase. Therefore this regime 
is particularly difficult to treat because on the one hand nucleonic degrees of freedom are no 
longer the proper degrees of freedom while on the other hand perturbative QCD calculations 
are notoriously difficult. 

For this reason we will evaluate the emissivity for unphysically high densities of n q ~ 
100 — 10 uq and extrapolate from this region down to the physical regime. At such high 
densities the momentum carried by the quarks is p ~ 1 — 40 GeV. This being in the per- 
turbative regime, we can estimate the strength of the QCD coupling constant a s employing 
the leading-log formula for the running coupling [derived from the /3-function to 0(«g)] 



O'.s 



El 

Air 



12tt 



(33-2iV F )m<g> 



(3.7) 



which gives a s ~ 0.48 — 0.13 — small enough for perturbative calculations. To check the 
effects of higher order correction we also use the formula 



Ait 



4tt 

(3 u 



2/3i1iim 401 



+ 



(31 u $hi 2 u 



i 1x2 
mu 

2, 



8(31 



(3.8) 
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derived from the /3-function to 4th order in a s . Here 

2 1Q S033 32^ (n 2 ) 

Po = U--N F , f3 1 = 51--N F , [3 2 = 2857 - —N F + — N 2 F , and M = ln^f, 

(3.9) 

where A rs 200 MeV is taken as the QCD breakdown scale and N F = 3 is the number of 
flavors. The formula in Eq. ( |3.8|) changes a s by about 18 — 37% compared to Eq. ( |3.7| ). 

It is instructive to compare the number of nucleons and quarks N q participating 
in scattering for the two cases of nucleon-nucleon scattering (for simplicity we assume here 
neutron matter) and qq scattering in the degenerate limit. For low temperature the baryon 
number density % is approximately proportional to pp. For nucleons % ~ Pf/(3vt 2 ), while 
for quarks n B « p 3 F /n 2 . The number of particles N N [N q ] participating in scattering is 
roughly proportional to the momentum-space volume of a shell of the Fermi sphere with 
radius p F and thickness dE « T leading to Nn « 8irp F dp [N q « 72np F dp] in the nucleon 
[quark] case. Making use of the energy- momentum relations pdp = M dE [pdp = p dE] this 
becomes N N « 7r(37r 5 ) 1/3 nB /3 MT [N q = 727r 7 / 3 n| /3 T] . Using M « 940 MeV the ratio of N N 
and N q for a density of n B = 10n is 



Nn 



0.3 (3.10) 



and hence, despite the fact that the quarks have a higher number of degrees of freedom, the 
number of particles participating in scattering is not much different in the two cases. 



A. Quark-Quark Scattering Cross Sections 



The matrix elements for the scattering of quarks to first order in a s can be calculated 
from the one-gluon exchange diagrams and are given in ||25|| P|. For quarks of different flavor 
the squared matrix element can be written as 

2 2 

\ A qq ,\ 2 = 256n 2 a f ^ (3.11) 
and for quarks of the same flavor as 

|^,|. = 26^04(^ + ^-13 (3 ' 12) 
(see Fig. ^). Here s, t, and u are the usual Mandelstam variables given in the COM frame 



3 These matrix elements, however, include an averaging [summing] over initial [final] spin and 
color states. Since for the calculation of the emissivity we need to sum over spin and color of both 
the initial and final states, we have to include an additional factor of 36 in the squared matrix 
elements. 
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FIG. 2. The lowest order diagrams for scattering of quarks of the same flavor. 



for zero quark mass by s = 4p 2 = 4p' 2 , t = —2p 2 (l — cos# cm ), and u = —2p 2 (l + cos# cm ), 
where 9 cm is the COM scattering angle. 

Note that similar diagrams exchanging a photon instead of a gluon are suppressed by a 
factor of ctQED / «s and need not be considered. 

The matrix elements in Eqs. (|3.11|) and ( |3.12 ) can become highly singular because of the 
t 2 and u 2 in the denominators. The momentum q transfered by the exchanged gluon vanishes 
as 9 cm — ■> or 9 cm — > it, which corresponds to colinear scattering. On the other hand, the 
energy-momentum tensor also vanishes for colinear scattering and one might think that 
this would render the singularity spurious. However, this is not the case. An analysis of the 
functions gx shows that for 9 cm [9 cm ~ n] gx behaves like ~ 9 cm 2 [~ (it — 9 cm ) 2 } while 
\A qq \ 2 and |^4 gg '| 2 behave like ~ l/9 cm 4 [~ — 9 cm ) 4 ]. 



B. Curing the Divergences 

In the case of electron-electron scattering in an electron gas the divergences that one 
encounters are of a similar origin and are usually taken care of by including effects which 
screen the bare particles arising from the interaction of the electrons and photons with the 
surrounding matter. We will employ, without justification, a similar procedure in the present 
case of a QCD plasma. 

The bare gluon propagator D^ u is related to the gluon self-energy by Dyson's equa- 
tion 



UJ 



q 2 ) + n 



(3.13) 



Following Ref. |]26|| , the matrix element for scattering of like quarks for small momentum 
transfer is equal to the one for unlike quarks and can be written as 



qq' 



128vr 2 a 2 



^1 — X 2 ) COS0 



q 2 + n 2 q 2 - uj 2 + lit 



(3.14) 



where II; and H t are the longitudinal and transverse gluon polarization functions given for 
q < »q by 



n, 



X . X + 1 

— In 

2 x - 1 



and 11/ 



x x(l — 



X ) . X 

m — 



x 



(3.15) 



Here x = u/q and qr> is the Debye wave number for cold quark matter of three flavors 
given by q 2 D = 6a s n 2 /7r 2 . Note that for quarks of equal mass in the COM frame u = and 
therefore U t = 0; the transverse part of the propagator is still unscreened, at least for U t 
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calculated to zeroth order in x. Since T <C fi the participating quark's momenta are near 
the Fermi surface and hence u ~ 0, independent of the chosen reference frame. The problem 
of an unscreened transverse gluon propagator can be solved in some cases by retaining 0{x) 
terms in Eq. (|3.15|) and is related to frequency-dependent screening |[27| . 

Because of these complications we adopt here a much simpler method to incorporate the 
complicated screening effects. We will render the amplitudes in Eqs. (|3.11| ) and (|3.12| ) finite 
by including a Debye screening mass mo, by substituting 



t -> t 



m 



D 



and u 



u — m 



D 



(3.16) 



for u and t in the denominator. For a small coupling constant the Debye mass to lowest 
order is given by 



ni 



D 



(3.17) 



which makes it about 60—20% of a typical momentum transfer for densities n q = 100 — 10 7 no- 
To consider higher order corrections to and to account for other many-body effects in 
the QCD plasma at these high densities we carry out the calculation for the emissivity 
not only using the value uid from Eq. (|3.17|) but also the values m£>/3 and 3m£>, thereby 
spanning about an order of magnitude in the cutoff mass. As an example we show in Fig. |3| 
the dimensionless function B{6 cm ) defined by 








30 



60 90 120 150 180 

9 Tdegreesl 

cm 

FIG. 3. The function B(6 cm ) for n q = 10 5 no and for different cutoff masses m£>/3 MeV (solid 
line), mo MeV (dashed line), and 3mD MeV (dot-dashed line) with uid = 2250 MeV. The typical 
COM momentum is 8400 MeV and a, is set to 1. 



\A 



1h 



+ 



\A, 



sin 8 r 



(3.18) 
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which enters the formula for the emissivity for a varying cutoff mass mp. The factor sin 2 9 CTa 
in B(6 cm ) is an approximation for Jq dxx d ~ l gx{x-, 1, # C m) as explained in the next section. 
B(6 cm ) has been calculated for a plasma at the density n q = 10 5 no where a typical COM 
momentum is about 8400 MeV and rs 2250 MeV^. The variation of mjj of about one 
order of magnitude causes a variation in B(8 cm ) of about 1-2 orders of magnitude. 

Other possible ways to regulate these divergences could be to introduce a cutoff of order 
mo for the gluon- momentum integration or to cut off the integration over 9 cra at some small 



angle 1/A when calculating observables. The latter technique was used by Weinberg [23 
to calculate the power produced in gravitational radiation from the hydrogen plasma in the 
sun. 



IV. GRAVITATIONAL EMISSIVITY OF A GAS OF RELATIVISTIC FERMIONS 

In order to calculate the emissivity due to graviton radiation from the two-body scattering 
reaction qq — > qqX in a gas of relativistic quarks we have to integrate over the phase space 
of incoming (p[, p" 2 ) and outgoing (p[, p' 2 ) quarks as well as that of the emitted gravitational 
particle (k). The formula for the emissivity is then given by 



dS 



x 



dt 



-j-j. d 3 Pi d z p\ 



.i=l. .2 



(2tt) 3 2£ 4 (2tt) 3 2^ 



/i/ 2 (l-/i)(l-^) 



A J 

x(2tt) 4 5V + p 2 -p[ -p' 2 - k)-£. (4.1) 

Here, again, the X stands for the type of gravitational radiation, massless 4D-gravitons 
(X = g), massive KK-gravitons(X = h), or massive KK-dilatons (X = <fi). The functions 
/ are Pauli blocking factors given by fi = (exp[(Ei — fi)/T] + furthermore = \pi\ 
for massless quarks. The summation over the KK-modes j doesn't apply to the massless 
4D-graviton, obviously. 

This formula applies for scattering of differently flavored quarks. A symmetry factor of 
1/4 would have to be included for scattering of quarks having the same flavor to account 
for identical particles in the initial and final states. We will consider this symmetry factor 
at the beginning of the next section. 

In the soft-radiation limit we can neglect the graviton momentum k in the momentum- 
conserving delta function. Furthermore, by introducing the total momentum P = p~[ + p~2 = 
Pi + P2 an d relative initial and final momenta p = (p~[ — p~2)/2 and p' = — p'2)/ 2 we 
can reduce the number of integrations by exploiting spherical symmetry and momentum 
conservation to get, after inserting Eq. (|2.21| ) for the massless quark case (y — 1), 



dFv f /*1 /*! r°° /*°° /*°° r, /*27T 

-w=w+vl*"L da » t L ia "'L dPpl l dppl l *' p I 



.1 



x/i/ 2 (l-/0(l-/D 



x + E 2 -E[- E' 2 )p 2 g x (x, 1, 6 cm ) \A(p, 6 cm )\ 2 . (4.2) 



EiE 2 E' 1 E' 2 
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Here we have also used the COM frame for the amplitude A. The angles 9 and 9' are defined 
by p ■ P = pP cos 9 and p' ■ P = p'P cos 9', while <fi is the angle between the projections of p 
and p 1 on the P-plane, that is 

cos 9 cm = cos <p sin 9 sin 9' + cos 9 cos 9'. (4.3) 

If the mass splitting of the KK-modes becomes comparable to the experimental energy 

— * 

resolution — which is true for about d < 6 — the summation over the modes j can be approx- 
imated by an integration: 

— ► R d & d f dmm d - 1 = R d Q d u d ^ dxx d -\ (4.4) 

J 

where d is the number of extra dimensions and fl d is the surface area of the rf-dimensional 
sphere given by tt d = 2-K d / 2 /T{d/2). 

To further simplify the integration we make another change of variables to the initial 
and final COM total energies E T and E' T given by 



P 2 P 2 
E T = E 1 + E 2 = d M 2 + — + p 2 + Ppcos9 + d M 2 + — + p 2 - Ppcos9 (4.5) 

and a similar expression for E' T . The Jacobian J of the variable transformation (p,p') — > 
(E T , E' T ) is given by 

dE T dE' T \ 2 ^ Et 2 _ 4M2 _ p 2(jB 2 _ p2 cog2 )3/2 J [ J 

(4.6) 

After interchanging the E T and i?^ integrations and shifting the integration variables the 
final result for KK-gravitons and KK-dilatons becomes 

dS x G 



dt QAtv 7 



x - N R d tt d / d9 / d& / dP da da' sin 9 sin 9'P 2 a d Jp 2 p' 2 
Jo Jo Jo Jo Jo 



x/i/ 2 (l - - ti) EiE P 2 e>e> 2 d( t>\A{PM 2 [dxx^gxix, l,9 cm ), (4.7) 



where we have defined 



a' = E' T - 2]j M 2 + ^- and a = = E T - E' T . (4.8) 

Note that only the function gx depends onx = m^/uj. Also, the integration over is hidden 
via 9 cm in the function gx and therefore nontrivial. 
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A. Approximation Formula for Degenerate Quark Matter 



Quark matter with a baryon number density ub ~ 10n has at a temperature T ~ 
50 MeV a quark chemical potential \i q ~ 500 MeV and is highly degenerate, \i q jT 3> 1. At 
higher densities matter is even more degenerate. Therefore one can assume that radiation 
only arises from scattering involving quarks near the Fermi surface in the initial and final 
states. Assuming soft radiation the energy of the quarks in Eq. (11 ) can be set to pp causing 
a decoupling of the energy and angular integrations: 



dS 



x 



dt 



^JdujJ dE x dE 2 dE[ dE' 2 hf 2 {l - f[)(l - f 2 )5{E x + E 2 -E[-E' 2 -lu) 

<Y { dtt t dtt 2 dn[ dtt'Siti + p 2 - p\ - p\) (4.9) 
. i aw 



The energy integral is well known [28, 29 1 and given by 

J dEt dE 2 dE[ dE' 2 fJ 2 (l - f[)(l - ti)6{E! + E 2 -E[-E' 2 -uj) 



u(u 2 + Aix 2 T 2 ) 



. (4.10) 



Replacing the sum over rrtj by an integration [see Eq. (4.4)1 we can carry out the ^-integration 
to get 



dS 



x 



Pf 



G N Q d R d Y d r da sin 3 a [* d6 cm \A\ 2 (9 cm ,p) [ dxx^gx&l, 9 cm ), (4.11) 
Jo Jo Jo 



dt 32tt 7 

where p = Pf sin a. The function Yd is given by 

rp4+d 



6 



47r 2 r(2 + d)C(2 + d) + r(4 + d)C(4 + d) 



(4.12) 



where ((x) is the Riemann Zeta function. This reduces the number of integrations to be 
carried out numerically to three. 

The temperature dependence of the emissivity in Eq. (|4.11|) can be easily understood. 
The degenerate quarks must have energies that lie within ~ T of the Fermi surface and 
are contributing one power of T each. Because of the energy conserving Delta function the 
energy of the graviton is of order T as well. Each extra dimension contributes one more 
power of T from the KK mode. Because only four out of the five energies are independent, 
the power of T gets reduced by one and we end up with the T 4+d dependence. 

Furthermore, for energy loss due to KK-graviton (X = h) radiation — the dominant 
mode — we can replace the x-integration by 



dxx d 1 gh(x, 1, 9 

cm ) 



(2 -2 log 2) sin 2 9 cm for d = 2 
0.38 sin 2 9 cm for d = 3 ' 



(4.13) 



which introduces an error of about 10% but reduces the number of numerical integrations 
to two. 
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V. RESULTS AND COMPARISON TO NUCLEON-NUCLEON SCATTERING 



For quark matter with equal number fractions for the three flavors u, d, and s, X u = 
Xd = X s = 1/3, the total emissivity £ can be written as a sum over possible combinations 
of two out of the three flavors in the plasma: 



£ — X! ^£qqXg + £ U( iX u X d + £ US X U X S + £ ds X d X s 

q=u,d,s 



(5.1) 



where the 1/4 is the symmetry factor for the scattering of like quarks. Then the formula for 



the emissivity [Eq. ( |4.11] )1 can be used to calculate the total emissivity by simply substituting 



Note that for massless quarks the trace over the energy-momentum tensor T^ v is zero and 
so the KK-dilatons decouple and need not be considered here. Furthermore, the contribution 
of the massless 4D-graviton can safely be neglected since it is several orders of magnitude 
smaller than that for the massive KK-gravitons. 

In Fig. |] we show the graviton emission rate from quark matter for a density range of 



10 



28 



de/dt 26 
[ergs/g/s] 10 




n [n Q ] 



FIG. 4. Emissivity from quark matter for high densities of 100 — 10 7 no- The solid [dashed] 
line shows the result using the formula for the QCD coupling constant derived from the /3-function 
to 0{a,g) [0{a s )]. The cutoff mass mn is calculated using Eq. fl3.17| ). The upper and lower 
dotted lines correspond to a cutoff mass of mjj/3 and 3m£>, respectively. The results for nuclear 



matter at lower densities from Ref. [16| are also included, using the approximations for degenerate 
(dot-dashed line) and non-degenerate matter (dot-dot-dashed line). All calculations are for the 
case of d = 2 extra dimensions and a temperature of T = 50 MeV; the size of the extra dimensions 
is taken to be 1 mm. The gray area denotes the region of realistic density of n q ~ 10 — 20no- 



100 — 10 7 no for the case of d = 2 extra dimensions and a temperature of T = 50 MeV, 
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assuming that the size of the extra dimensions is R = 1 mm. The solid line shows the 



result using the cutoff mass calculated from Eq. ( 3.17 ) valid for small a s and using 
the formula Eq. ( |3.8| ) for the QCD coupling constant derived from the /3-function to 
To estimate the error introduced by calculating only first order diagrams we give also the 
result for the emissivity using the (larger) a s [see Eq. ( |3.7| )1 calculated from the /3-function 
to 0(a 2 s ) (dashed line). This results in an increase in the cross section which gets partially 
compensated by a larger cutoff mass. Therefore the emissivities for both values of a s are 
very similar. 

To account for uncertainties in the cutoff mass (higher order corrections, unaccounted 
many-body effects) we also show the emissivity using (upper dotted line) and 

3m d (lower dotted line). With mp spanning about one order of magnitude in the cutoff 
mass the emissivity varies by about 1-2 orders of magnitude. We also present the emission 



rates from nucleon-nucleon scattering |16| for lower densities. 

Note that it is not sensible to compare the emission rates for quark matter and nuclear 
matter for an equal density, since for densities much higher than 5no nucleons are certainly no 
longer the appropriate degrees of freedom. Similarly, for densities smaller than about lOOno 
the QCD coupling constant becomes too large so that perturbative QCD calculations will 
no longer be applicable. It is therefore necessary to extrapolate from the two regions where 
reliable calculations are possible to the interesting transition region with n g = 10 — 20n 
(gray area in Fig. |4]) . 

Determining the emissivity for the physical density region from Fig. f| for d = 2 and 
carrying out a similar analysis for d = 3, we can now use the Raffelt criterion, limiting the 
energy lost into any non-standard physics channel to 10 19 ergs/g/s, to calculate bounds on 
the size of new extra dimensions. The radius R enters the calculation through the prefactors 



in Eqs. ( |4.7| ) and ( flllp . Our bounds are calculated to be 



R < 2.9 x 10~ 4 mm for n = 2 and (5.3) 
R < 3.9 x l(r 7 mm for n = 3 (5.4) 



which are quite similar to the ones determined for nuclear matter |L6|. 

It should be kept in mind that these bounds become somewhat weaker [stronger] as 
tub is taken to larger [smaller] values thereby decreasing [increasing] the emissivity. Our 
conservative estimates are therefore weaker than those found from previous nuclear matter 
calculations. 



VI. CONCLUSION 

Supernovas like SN1987a provide one of the most stringent constraints on the size of 
gravity-only extra dimensions. Bounds on these sizes have been calculated before [I|-T7 



assuming a density of ~ 1 — 3n and a typical temperature of T « 50 MeV where matter is 
nucleonic. The relevant degrees of freedom for this regime are protons and neutrons. 

In view of the uncertainties with which the parameters governing the condition of matter 
at the star's inner core, such as density and temperature, are given, and in view of the fact 
that, as can be found from various nuclear matter equations of state, matter in the star's 
core lies very close to a phase transition to either a deconfined QCD plasma or, for lower 
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temperatures, a color superconducting phase, it is worthwhile and important to investigate 
how emission rates and bounds on the radius of the extra dimensions change in these regimes. 

In this paper we have calculated emissivities and bounds on the size of gravity-only extra 
dimensions from a deconfmed quark-gluon plasma in the core of a star seconds after the 
supernova core bounce. We have not considered the case where the temperature falls below 
the critical temperature T c defining the transition of quark matter to a color superconducting 
state. 

By extrapolating the emissivity from a quark-gluon plasma at very high density (n q = 
100n — 10 7 n Q ) down to the physically interesting region of n q ~ 10 — 20n we find that the 
KK-graviton emissivity is comparable to those found from nuclear matter calculations. We 
have examined the error introduced into our calculation due to the fact that we are working 
only to first order in the QCD coupling constant. We estimate this error to be at the most 
100%. To consider the rather uncertain value of the cutoff mass we have calculated the 
emissivity using a cutoff mass which varies about one order of magnitude about the central 
value from the lowest order perturbative calculation [Eq. (ft.l7|) l- We find that this causes 
an uncertainty in the emissivity of about two orders of magnitude. 

We have not examined any other many-body effects which could modify the vacuum 
rates for KK-graviton emission, such as the Landau- Pomeranchuk-Migdal effect, which would 
decrease the emissivity because of multiple quark scattering before the graviton gets emitted. 

Note that the source of uncertainties is not the soft radiation theorem, which is well 
suited for the degenerate case, but rather our lack of knowledge about the nuclear physics 
nature of QCD. This deficit of understanding to only impacts our calculation of the quark- 
quark scattering cross section and the treatment of the QCD plasma, but also forces us to 
access the interesting transition region only with the aid of extrapolation. 

Despite the significant uncertainties in our calculation for quark matter, we believe to 
have ruled out the possibility of a significantly larger emissivity compared to the nuclear 
matter case. It is up to further investigation to provide better information on the conditions 
which govern the state of the matter at the star's core during and shortly after a supernova 
event. It is unlikely that the whole core will consist of dense quark matter, more likely is 
a scenario where at the core's center matter is deconfmed and undergoes a phase transition 
to nuclear matter as one goes towards the surface. Furthermore, a better understanding of 
QCD in the non-perturbative regime and at finite density would enable a reduction of the 
uncertainties. 
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APPENDIX A: RESULTS FOR THE ANGULAR INTEGRATION 



The functions gx(x,y,9 cm ) defined in Eq. (|2.2Up are calculated for the three cases of 
4D-graviton (g), KK-graviton (h), and KK-dilaton ((f)) to be 

1/ 1-L 2 , (l + 6y 2 + y 4 )log^i 2y 2 + y* + 4y 2 cos9 cm + cos29 cm 
9g = 7\- 1 +y 



4[ 2y W(x,y,cos0 cm ) 

V - l) 2 + 2(y 2 + cos 6 cm ) 2 + 2{y 2 + cos cm ) W(x, y, cos 9 C 



x log 



(ya - If 



+7<cos6» cm -> -cos0 cm l, (Al) 



1 

6 



0h = -Vl^.r 



4(, 2 -l) 2 , 2(^ + V + l)log^_ 



x 2 + y 2 - 1 vi - x 2 y 

12y 2 cos # cm + (2y 2 + l) 2 + 3 cos 26 cm 



(x 2 — l)W(x, y, cos 9 C 



x log 



{x 2 - 1) 



(x 2 +y 2 -l)' 2 
(x 2 -l) 



+ 2 + cos # cm ) 2 + 2 + cos 9 cm ) W(x, y, cos 9 C 



(x 2 + y 2 - 1) : 



+ - Vl — x 2 \ cos # c 
6 



COS 61. 



(A2) 



and 



d- 1 
3(rf + 2) 



Vl - x 2 (y 2 - 1)< 



a; 2 + y 2 — 1 



1 (2/ 2 -l)lo. 
+ - 



y+Vl—x' 2 
y—\/l—x 2 



2V1 - ^ 2 Z/ 



+ 



^s 2 — 1) W(a;, cos # c 



x lot 



(x 2 - if 



(x 2 +y 2 -l) 2 



(x 2 -l) 2 



+ 2 + cos 9 cm ) 2 + 2 + cos # cm ) W(x, y, cos 9 C 



{x 2 + y 2 -l) 2 



d-1 



Vl — x 2 (y 2 — 1) < cos 6* cm — > — cos 9 C 



3(d + 2) 

The function W(x, y, h), introduced for convenience, is defined as 



(A3) 



W(x,y,h) 



\h + l)[-a; 2 - 2y 2 + (x 2 - l)h + 1] 
x 2 — 1 



(A4) 
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